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Abstract 



We search for a new light non-Standard Model CP-odd Higgs boson, A , decaying to 
t + t~, and invisible in radiative decays of the T(2S) and T(3S). We search for the 

light scalar particle produced in single-photon decays of the T(3S) resonance through the 
process T(3S) — > 7 A , A — > invisible. We also search for evidence of an undetectable T(1S) 
decay recoiling against the dipion system. The data samples contain 99 x 10 6 Y(2S) and 
122 x 10 6 T(35) decays collected by the BABAR detector at the SLAC PEP-II B factory. 
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1 Introduction 

The search for the origin of mass of fundamental particle is a cornerstone of the Standard Model 
(SM). The Higgs mechanism is one of the great interests in the particle physics. In recent years 
many motivations for searching a Higgs boson. Direct searches at the CERN e + — e~ Collider 
(LEP Experiment) via the process e + — e~ — > ZH for the Standard Model Higgs boson, a CP- 
even scalar, has excluded a lower bound on its mass below 114.4 GeV pQ. The search for the 
SM Higgs particle is continued at the Tevatron (The CDF and DO Collaborations) for a Higgs 
boson decaying to W + W~ and resulting limits on Higgs boson production exclude a SM Higgs 
boson in the mass range 162-166 GeV at the 95% C.L [2J. The search is currently continued at 
Large Hadron Collider (LHC) with the relative increase in charged-particle multiplicity from 
^ = 0.9 to 7 TeV [3S]. 

In recent years a number of theoretical models predicted for the existence of a light CP-odd 
Higgs boson (A ) related to the Next-to-Minimal Super symmetric Model (NMSSM). Direct 
searches constrain the mass of m^o < 2rrib, where is the b quark and the decay of A — > bb 
is forbidden [5]. Of particular interest is to search the lightest CP-odd Higgs boson in Upsilon 
decays such as T(1S) and T(3S). In this paper we present a search of the CP-odd Higgs boson 
decaying to t + t~, fi + fj>~, and invisible. The large data sets available at BABAR experiment 
allow us to place stringent constraints on such theoretical models. 
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2 Data Sample 



The data used in this paper were collected with the BABAR detector at the PEP-II asymmetric- 
energy e + — e~ storage ring at the SLAC National Accelerator Laboratory , operating at T(2S) 
and Y(3S). For searching for A decays to either r + r _ or fi + fi~, we use a data sample of 
122 x 10 6 T(3S) events that corresponds to an integrated luminosity of 28 fb . We also used 
a data sample of 79 fb -1 accumulated on the T(AS) resonance (ON-resonance) and 8 /6 _1 49 
MeV below the T(4S) resonance (OFF-resonance) for studying the background. 

For searching an undetectable T(1S) decay recoiling against the dipion system we used 
the data were taken using an upgraded muon system, instrumented with both resistive plate 
chambers and the limited streamer tubes between between steel absorbers. The data triggers 
was modified to substantially increase the pion trigger efficiency. The data sample of 96.5 x 1CP 6 
T(3S) were used. The BABAR detector is described in detail elsewhere [6]. 

3 Analysis Method 

3.1 Search for A in the decays of Y(3S) ->- -yA°, A t+t~ 

We search the light CP-odd Higgs boson, ^4°, via the decays T(3S) — > jt + t~ for a wider mass 
range, 4.03 < m T+T - < 10.10 GeV/c 2 [7], over the mass ranges by the CLEO Collaboration [8] 
and the DO Collaboration [9], respectively. We exclude the mass region of 9.52 < m T + T - < 9.61 
GeV/c 2 because of the irreducible background photons produced in the decays of T(3S) — > 
7Xw(2P), Xbj(2P) ^(lS), where J = 0,1,2. We scan for peaks in the distribution of the 
photon energy, P 7 , corresponding to peaks in the rr invariant mass as given by 

m l+r- = m T(3S) ~ 2m r(3S)Ey, (1) 

where m T ^ s ^ is the T(3S) mass (10.355 GeV/c 2 ) and the P 7 is measured in the T(3S) center- 
of-mass frame. 

We select events in which both r-leptons decay leptonically to either r + — > e + v e v T or r + — > 
[i + v^v T . The events are then required exactly two charged tracks to reduce the background and 
to contain at least one photon with energy > 100 MeV in the electromagnetic calorimeter. In 
addition both charged tracks are required to be identified as leptons either electron or muon. 
The residual background is mostly due to radiative decays of e + e~ —> ^t + t^ and higher order 
of QED processes, including two-photon reactions such as e + e~ — > e + e~e + e~ and e + e _ — > 
e + e~fi + fi~. However this residual background has smaller contributions compared to other 
T(3S) decays and the continuum background from the non-resonance T(3S) decays of e + e~ — > 
7* — > qq, where q = u,d, s, c. 

The events T(3S) — > 7Xw(2P), xw(2P) — > jY(nS), and T{nS) decays to t + t~, where 
J = 0,1,2 and n = 1, 2 are expected to peak in the photon energy distribution when the 
photon comes from the decays of x^j(2P) — > jT(nS) is misidentified as the radiative photon 
from the T(3S) decays. We use a Crystal Ball function (CB) [10] to describes each of the peaks 
that comes from to the decays of Xfoj(2P) — > jT(1S) in the photon spectrum. The mean values 
for the xw(2P) and the Xfei(2P) CB function are fixed to the PDG values [TT], and the width 
values are fixed to the Monte Carlo (MC) resolution, but the mean and the width of X62(2P) are 
free parameters. The search for the signal Y(3S) —¥ jA°,A° — > t + t~ is performed by scanning 
for peaks in the photon energy distributions. Figure [T] shows the fits to the photon energy, P 7 , 
distributions in the different t + t~ decay modes fitted in the region 0.2 < E gamma < 2.0 GeV. 

No evidence for a signature of the light Higgs boson decays to r pairs is observed. We 
measure the branching fraction of B(T(3S) jA°) x B(A° t+t") < (1.5 - 16) x 10~ 5 90% 
confidence level (C.L.). We also set a 90% C.L. upper limit on the t + t~ decays of the rjb at 
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Figure 1: The photon energy distributions for the different rr-decay modes. The data is shown 
in filled circles, the dotted lines represent the contributions from T(3S) — > , yxbj(2P),Xbj(2P) — > 
7T(25 r ), dotted-dashed lines show the contributions from T(35) — )• 7x&j(2P), X6j(2P) — >■ 
7T(1S'), and the solid lines show the total background functions. The different between the 
background function and the data divided by the uncertainty in the data for each rr-decay 
modes is shown in b, c, and d. 



3.2 Search for A in the decays of T(3S) ->■ 7 A , A fi+fj,- 

We search for a resonance in the ji + [a~ invariant mass distribution for the fully reconstructed 
final state T(2S, 35) — > 7vl ,yl — > fi + fi>~. We assume the ^4° resonance is a scalar (or pseu- 
doscalar) particle and the width of the resonance is negligibly small compared with the experi- 
mental resolution where the mass of light Higgs boson, m^o, is sufficiently far from the mass of 
the Vb H2J. 

We select events with exactly two oppositely charged tracks and a single energetic photon 
with a center-of-mass energy E* < 0.2 GeV, while allowing additional photons with the center- 
of-mass energies below 0.2 GeV. We assign a muon mass hypothesis to the two tracks and 
require that at least one is positively identified as a muon. We then perform a kinematic fit to 
the T candidate from the two muon candidates and the energetic photon. 

The backgrounds are dominated by the QED processes: the continuum background e + e~ — > 
7^ + /U _ and the initial-state- radiation (ISR) production of p°, <p, J/ip, ip(2S), and T(1S) vector 
mesons. We suppress the background contributions from the decays of p° — > 7r + 7r~ final state 
in which the pion is misidentified as a muon by requiring both charged tracks are positively 
identified as muons in the range 0.5 < m^o < 1.05 GeV. We also suppress the decays of 
T(2S) -> 72X&(1P), XbO-P) -> 7iT(lS) [^(35) -> 72X6(2P), Xb(2P) 7iT(lS)}, where 72 is 
the secondary photon by requiring that no secondary photon above a center-of mass energy of 
0.1 GeV [0.008 GeV] is present in the event. 
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Since we do not observe a significant excess of events in the range of 0.212 < m^o < 9.3 
GeV, we set upper limits on the branching fractions of B(T(2S) — > jA ) X B^ and B(T(3S) —> 
jA°) xBuu, where Baa is the branching fraction of A — > \x\i. The limit vary from 0.26— 8.3 x 10~ 6 
for B 2S and 0.27 - 5.5 x 10" 6 for B 3S - The 90% C.L. Bayesian upper limits are computed with 
a uniform prior and assuming a Gaussian likelihood function. We set an upper limit on the 
coupling f${m A o = 0.214 GeV) < 1.6 x 10~ 6 at 90% C.L., by assuming W = 1. This result 
is significantly smaller than the value required to explain the HyperCP events as light Higgs 
production |13[ I14j. Figure [5] shows the upper limits on the branching fractions as a function 
of the mass m^o . 
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Figure 2: The 90% C.L. upper limits on the branching fractions of (a) B(T(2S) — > 7^4°) x Baa, 
(b) B(T(3S) — > jA°) x Bpn, and (c) effective coupling fj- x £> MM as a function of m A o (color 
online). The shaded ares show the regions around the J/ijj and tp(2S) resonances excluded from 
the search. 



We also measure the branching fractions of B(T(2S) — > 7775) x B(rjb — > fi + fi ) = (—0.4 ± 
3.9 ± 1.4) x 10~ 6 and B(T(3S) -> 7%) x B(rj b -> n+n~) = (-1.5 ± 2.9 ± 1.6) x 10~ 6 , where the 
first uncertainty is statistical and the second is systematic, dominated by the the uncertainty in 
r(%)- Using the BABAR measurement of B(T(2S) — > 7%) and B(T(3S) — > 7%) [12], we derive 
B(r) b -> /i + /x~) = (-0.25 ± 0.51 ± 0.33)% and B(t] b -> /i+ju") < 0.9% at 90% C.L. This limit is 
consistent with the mesonic interpretation of the rjb state [15J . 

3.3 Search for Y(3S) 7 A , A Invisible 

We search for the light scalar particle produced in single-photon decays of the T(3S) resonance 
through the process T(3S) — > "fA°, A — > invisible. We split the data sets into two broads energy 
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ranges based on the energy of the highest-energy in the center-off mass in our electromagnetic 
calorimeter cluster. The high-energy region is 3.2 < E* < 5.5 GeV, where E* = E com /2. The 
backgrounds in the high-energy region are dominated by the QED process such as e + e~ — > 
77, where the photon energy distribution for this process peaks. The low-energy region is 
2.2 < E* < 3.7. The backgrounds in this region are dominated by the radiative Bhabha events 
e + e~ — > je + e~ — > in which both electron and positron are misidentified. 

We extract the yield of signal events using a likelihood fits to the distribution of the observ- 
able missing mass squared in the mass interval < m A o < 6 GeV: 

m\ = 7t4 (3S) - 2E*m T(3S) , (2) 

where mx is the missing mass squared, mj-(3S) is the T(3S) mass, and E* is the photon energy 
in the T center-off-mass system. The E* is given by 

The current precise limit on the measurement of the branching fraction of B(T — > jX) with 
X — > invisible is given by CLEO collaboration on T(1S) [16j . 

In the high-energy range, the selection efficiency for signal is 10 — 11% depending on the mass 
m A o, and it is below 10 -5 for e + e~ — > 77 events, whereas in the low-energy range the selection 
efficiency for signal is about 20%. The signal PDF is described by the Crystal Ball function [10] 
that centered around the expected value of m? x = m^A . We compare the distributions of the 
simulated and reconstructed e + e~ — > 77 events to determine the the uncertainty of the PDF 
parameters. 

We determine the PDF as a function of mass m A o. In the low-statistics simulated samples 
of signal events, we exclude 6 < m A o < 7.8 GeV and in the high-statistics simulated samples of 
signal events, we excluding the values of m A o < 6 GeV due to low statistics. The fit results give 
us signal yields of N s i g = 119 ± 71 (1.7a) and N S i g = 37 ± 15 (2.6a) for the low-energy dataset 
and the high-energy dataset, respectively. The fit results on the low-energy range (m^o = 7.275 
GeV) and the high-energy (m^o = 5.2 GeV) are shown in Fig. [3l 

We do not observe a significant excess of events above the background in the mass range of 
< m^o < 7.8 GeV., and we set 90% C.L. upper limit on the branching fraction of B(T(3S) — > 
j A ) x B(A° -> invisible) at (0.7 - 31) x 10~ 6 in the mass range of m A o < 7.8 GeV [IT]. These 
results are preliminary. 



3.4 Search for T(1S) Invisible 

The nature of dark matter is one of the most challenging issues in the elementary particle 
physics. Observation of the SM particles coupling to undetectable (invisible) particles might 
provide information on dark matter candidates. The dark matter couples to the SM through 
the mediating boson. This particle is expected to be a vector, scalar, or pseudoscalar. The dark 
matter to show up dominantly in interactions with heavy fermions such as b quarks \18\ [T9] . 

We select events in the invisible subsample by requiring exactly two tracks originating from 
the interaction point (IP) with opposite charge. These track have a center-off-mass momentum 
of p* < 0.8 GeV/c and they are consistent with pions from dipion transition with the invariant 
mass of 0.25 < < 0.95 GeV/c 2 . The dipion recoil mass is given by 

Mi c = s + Ml-2 v r {s)E * n) (4 ) 

where M rec is the reconstructed dipion recoil mass, E* n is the center-off-mass energy of the 
dipion system, and y^(s) = 10.3552 GeV/c 2 . The overall efficiency of the selection events is 
about 64%. 
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Figure 3: [Left] Sample fit to the low-energy dataset (83 x 10 6 T(3S)) decays. The bottom plot 
shows the data (solid points) overlaid by the full PDF curve (solid blue line), signal contribution 
with m A o = 7.275 GeV (solid red line). [Right] Sample fit the high-energy dataset (122 x 
10 6 T(3S)) decays. The bottom plot shows the data (solid points) overlaid by the full PDF 
curve (solid blue line), signal contribution with m^o = 5.2 GeV, e + e~ — > 77 (dot-dashed green 
line), and the continuum background (black dashed line). 



We collect the invisible subsample by selecting tree-track and four-track events using the 
dipion system. We search high-momentum tracks from T(1S) decays by requiring only one or 
two additional tracks originating from the IP point, each with momentum p* > 2.0 GeV/c 2 . We 
then use our lepton identification to treat these tracks. If either both tracks passing electron- 
identification criteria the both tracks are treated as electron candidates; otherwise, both are 
treated as muon candidates. Figured] shows the M rec distribution from our invisible subsample. 

We determine the parameters of the PDF by fitting the M rec distribution in the four-track 
data subsample. We extract the peaking contribution by using an extended unbinned maximum 
likelihood fit and the non-peaking contribution using a first-order polynomial function. The 
peaking background is estimated using Monte Carlo subsamples of 1019 T(1S) — > e + e~ events, 
1007 T(1S) events, 92 T(1S) t+t" events, and 2.9 ± 1.3 T(1S) hadrons events. 

We search for evidence of the decays of T(1S) into undetectable final states recoiling against 
the dipion system in the T(3S) — > -k + it~Y{IS), using a sample of 9.14 x 10 6 T(3S) mesons. We 
find no evidence for the decays of T(1S) —> invisible. We set an upper limit on the branching 
fraction of B(T(1S) invisible) < 3.0 x 10~ 4 at the 90% C.L. [20]. 

4 Summary 

We search for a new light non-Standard Model CP-odd Higgs boson in the T(35), T(2S), 
and T(1S) systems. We also search for invisible decays of the T(1S) system. BABAR recent 
measurements are important input to rule out much of the parameter space allowed by the 
light non-Standard Model CP-odd Higgs and other models. These measurements also provide 
new and stringent constraints on the light CP-odd Higgs boson. In addition the measurement 
of the branching fraction of T(1S) —> invisible is extremely important as a sensitive probe of 
dark matter. 
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Figure 4: The maximum likelihood fit to the dipion recoil mass for data in the invisible sub- 
sample. Dash-dotted line is the nonpeaking background and the solid fitted-line is the peaking 
background. The solid line is the total fit function. 
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